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Tagged particle in a sheared suspension: effective temperature 
determines density distribution in a slowly varying external po- 
tential beyond linear response 
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Abstract. - We consider a sheared colloidal suspension under the influence of an external po- 
tential that varies slowly in space in the plane perpendicular to the flow and acts on one selected 
(tagged) particle of the suspension. Using a Chapman-Enskog type expansion we derive a steady 
state equation for the tagged particle density distribution. We show that for potentials varying 
along one direction only, the tagged particle distribution is the same as the equilibrium distri- 
bution with the temperature equal to the effective temperature obtained from the violation of 
the Einstein relation between the self-diffusion and tagged particle mobility coefficients. We thus 
prove the usefulness of this effective temperature for the description of the tagged particle behav- 
ior beyond the realm of linear response. We illustrate our theoretical predictions with Brownian 
dynamics computer simulations. 



Introduction. — While the principles of equilibrium 
statistical mechanics are well established, a general frame- 
work of its non-equilibrium counterpart is still lacking. 
Recently, there has been a lot of interest in various facets 
of non-equilibrium statistical mechanics [1]. One of the 
recurring concepts is the so-called effective temperature 
[5] . Originally, it was defined through the violation of the 
fluctuation-dissipation relation in mean-field spin systems 
[3] . Subsequently, a variety of definitions of the effective 
temperature, in terms of violations of different dynamic 
and static linear response type relations, have been inves- 
tigated using computer simulations of model fluids [2ll^{7j . 
Notably, there have been few theoretical studies of the ef- 
fective temperature [8UTT]. 



The original definition of the effective temperature and 
almost all subsequent investigations mentioned above were 
concerned with the linear response regime. On the other 
hand, in equilibrium statistical mechanics the temperature 
plays an important role also outside the realm of linear re- 
sponse. Hayashi and Sasa [8] were the first to discuss an 
application of the effective temperature outside the linear 
response regime. They studied an exactly solvable model, 
a single Brownian particle in a tilted washboard poten- 
tial. They showed that the effective temperature defined 



through the violation of the Einstein relation determines 
the density distribution in a long-wavelength external po- 
tential. Subsequently, there were three simulational stud- 
ies of the effective temperature outside the linear response 
regime. Importantly, these studies dealt with strongly in- 
teracting many-body systems. Zamponi et al. |I12| showed 
that the effective temperature enters into the fiuctuation 
relation jl3j for a driven glassy fluid. Ilg and Barrat jTl] 
showed that the effective temperature determines the bar- 
rier crossing rate in a driven glassy system. Haxton and 
Liu [15] investigated the usefulness of the effective tem- 
perature for the description of material properties (specif- 
ically, the shear stress and the average inherent structure 
energy) of driven glassy systems. To date, there is no the- 
oretical analysis of the effective temperature of a many- 
body system, outside the linear response regime. 

Establishing the usefulness of the effective temperature 
outside the linear response regime is essential because, in 
principle, the violation of the fluctuation-dissipation re- 
lation does not have to be interpreted in terms of the 
effective temperature. The best illustration of this fact 
is the difference between two studies of the same exactly 
solvable model, a single Brownian particle in a tilted wash- 
board potential. Hayashi and Sasa [Sj used the violation of 
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the fluctuation-dissipation relation to define the effective 
temperature. On the other hand, Speck and Seifert [TU] 
showed that the fluctuation-dissipation relation can be ef- 
fectively restored if fluctuations of the Brownian particle 
velocity are measured with respect to the local mean ve- 
locity. We should also mention a recent investigation [TB] 
which strives to provide a unifying framework for non- 
equilibrium fluctuation-dissipation relations. It is simi- 
lar to Speck and Seifert 's approach in that it focuses on 
an additive correction responsible for the violation of the 
fluctuation-dissipation relation. It de-emphasizes the ef- 
fective temperature altogether [T7]. These examples show 
the need for additional, fundamental understanding of the 
meaning of the effective temperature, especially beyond 
the realm of linear response. 

In this Letter we show theoretically that for an inter- 
acting many-body system one specific effective tempera- 
ture plays the role of the usual temperature outside the 
linear response regime. We consider a sheared colloidal 
suspension. In such a suspension the Einstein relation be- 
tween the self-diffusion coefficient and the tagged particle 
mobility is violated |18j and this violation can be used 
to define the effective temperature [HinillT]. We assume 
that the suspension is under the influence of a spatially 
varying external potential which changes slowly along one 
direction in the plane perpendicular to the flow and acts 
only on one selected (tagged) particle of the suspension. 
We show that the functional form of the tagged particle 
density distribution is then the same as that of the equilib- 
rium distribution. However, in this distribution the tem- 
perature is replaced by the effective temperature obtained 
from the violation of the Einstein relation along the same 
direction. Consequently, we can also prove that the same 
effective temperature is obtained from the tagged particle 
static linear response relation and the Einstein relation 
violation. 

Since the derivation of the main result is a little tech- 
nical, we flrst present the results and illustrate them with 
Brownian dynamics computer simulations. We sketch the 
derivation in the latter part of this Letter. 

Tagged particle density distribution. The A^- 

particle steady state probability distribution P^ describ- 
ing a suspension undergoing a shear flow and under an 
influence of an external potential satisfles the following 
equation, 

[n, + ar, • fio (9r,l^°'^*(ri))] P/(ri, ...,r^) = 0, (1) 

where Xlg is the Smoluchowski operator with shear flow. 



f^s = 



N 
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i=l 


V 


»#i=i 



F(r,,)-v(r, 



(2) 



In the above equations, fiQ and Do is the mobility and the 
diffusion coefflcient, respectively, of an isolated Brownian 
particle, F{rij) is the force acting on particle i due to 



particle j, r^ — Vi — Vj, and v is the uniform shear flow, 
v(ri) = r • r.i, with T being the velocity gradient tensor. 
Finally, V°^^ is an external potential. It acts only on the 
tagged particle and varies in the plane perpendicular to 
the flow, Sri ■ry<='=*(ri) = 0. 

We should emphasize that Eqs. ([TJ^l) without the exter- 
nal potential have been widely used in investigations of the 
dynamics and rheology of sheared suspensions [12] . There 
are two assumptions implicit in writing down these equa- 
tions. First, hydrodynamic interactions are neglected. 
Second, it is assumed that the shear rate is large enough 
to drive Brownian particles out of equilibrium, but it is 
small enough so that the solvent in which these particles 
are suspended is in equilibrium. The latter assumption 
implies that the Einstein relation for an isolated particle 
is satisfled, Dq — k^TiiQ, where T is the solvent temper- 
ature. 

In principle, in order to find the steady state tagged 
particle distribution rig we need to solve Eq. ([T]) and then 
integrate over the positions of all other particles. 



ns(ri) = / dr2...drArPg^(ri,...,rjv). 



(3) 



In the absence of the shear flow, we know the solution of 
Eq. ([T]) and we can easily flnd the equilibrium distribution 
for the tagged particle in external potential V''^^, 



n(ri)ocexp(-y<=''*(ri)/(A;Br)). 



(4) 



In the latter part of this Letter we use a Chapman-Enskog 
type expansion to show that, in the presence of the shear 
flow, for an external potential which varies slowly in the 
plane perpendicular to the flow, the steady state tagged 
particle distribution satisfies the following equation, 

9r, -{B-dr^+fi- (dr^V'^'ir^))) n,{r^) = 0. (5) 

Here rj^ is a two-dimensional vector in the plane perpen- 
dicular to the flow, and D and fi are the self-diffusion and 
the tagged particle mobility tensors which depend on the 
suspension's density, temperature, and shear rate. 

It follows from Eq. ([5]) that, if external potential y^''* 
varies only in the direction of unit vector n (in the plane 
perpendicular to the flow) , the tagged particle distribution 
has the form of the equilibrium distribution (|4]), 



(rj. • n) oc exp(-y°'^'(ri • h)/ikBT?")) 



(6) 



but with the temperature being replaced by the effective 
temperature defined in terms of the violation of the Ein- 
stein relation along direction n, 



T- 



■off 



n • D n 
h ■ fi ■ h 



(7) 



For large enough shear rates the system can become struc- 
turally anisotropic even though the density is still uniform. 
In particular, structure factors along different directions 
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in the reciprocal space vary differently with the shear rate 
[7]. In principle, one cannot expect that in such cases the 
effective temperature is isotropic. On the other hand, for 
a system with very slow dynamics {e.g. for a glassy sys- 
tem) there is a range of shear rates for which the system is 
strongly out of equilibrium but still structurally isotropic, 
with essentially shear-rate-independent structure factors 
along all directions in the reciprocal space [4]. In such 
cases it has been found that the effective temperature is 
isotropic in the linear response regime [4j. We expect (al- 
though we did not investigate it) that the same will be true 
also beyond the realm of the linear response. Of course, 
for any ergodic system T°^ — 5- T in the limit 7 — > 0. 

One of the consequences of Eq. ([6]) is the following re- 
lation between an infinitesimally small external potential 
^yext g^j^^ ^j^g resulting change Jng of the tagged particle 
distribution, 



Jns(r_L -n) 



kuT2 



SV^'^r^-n) 



(8) 



Eq. ([8|) is a generalization of the familiar equilibrium lin- 
ear response relation. It is analogous to static linear re- 
sponse type relations that were used to define effective 
temperatures by A. Liu and collaborators [S]. Impor- 
tantly, in Eq. (|8]) effective temperature (O appears in 
the place of a static linear response-based effective tem- 
perature. This is the first time that the effective tem- 
perature defined in terms of the violation of a dynamic 
fluctuation-dissipation relation has been shown theoreti- 
cally to be identical to the effective temperature defined 
in terms of the violation of a static linear response relation. 

Brownian dynamics simulations. We performed 
Brownian dynamics computer simulations as described in 
Ref. [7]. The system consists of A^ = 1372 particles inter- 
acting via a screened Coulomb potential, 



V{r) — Acxp (— K(r — a)) /r. 



(9) 



with A = Al^kBTa and kg = 24. We simulated the sys- 
tem at a dimensionless density, Na^/L^, equal to 0.408 (L 
is the box length), which corresponds to the hard-sphere 
volume fraction equal to 0.43 |5D|. We applied shear flow 
in the x direction with the velocity gradient in the y direc- 
tion, r — 7xy with 7 being the shear rate and a the unit 
vector along the a axis. In the following we use reduced 
units: we measure lengths in units of a, time in units of 
a^/Df), energy in units of fc^T, and effective temperature 
in units of T. 

To illustrate the validity of Eq. ([6]) we introduced a 
slowly varying external potential acting on particle num- 
ber 1 [2l]. We chose the simplest, plane wave potential. 



F'=''t(ai) = K)Sin(27rai/L), 



(10) 



where a — y,z. Note that 2Tr/L is the smallest wave- 
vector allowed by the periodic boundary conditions. We 
used Vo = 2, which is well outside the linear response 
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Fig. 1: (Color online) Tagged particle steady state density dis- 
tribution for suspensions with the external potential varying 
along (a) velocity gradient direction and (b) vorticity direction. 
Symbols: Brownian dynamics simulations' results. Solid lines: 
fits to a normalized version of Eq. (|6]). Dashed line: tagged 
particle density in the absence of the external potential. 



regime (see Fig. [T|). For each choice of the potential we 
performed 4 independent runs at the following shear rates: 
0, 5, 10, and 20. We monitored the tagged particle density 
distribution, ns{a), a = y,z. All results presented are 
averaged over 4 runs. 

In Fig. [1] we show tagged particle steady state density 
distributions. Qualitatively, with increasing shear rate 
the influence of the external potential decreases, which 
could be interpreted as an increasing effective tempera- 
ture. More importantly, the simulations' results can be fit- 
ted very well using equilibrium-like distributions ([6]) with 
an effective temperature being the only fit parameter. In 
Fig. [2] we show that the values of T?,*, a. — y,z, obtained 
from the fits agree very well with the previously obtained 
effective temperatures defined through the violation of the 
Einstein relation [7 . 

At the largest shear rates we found previously that the 
system is structurally anisotropic |7]. Thus, in principle, 
different effective temperatures can be obtained along dif- 
ferent directions. We note that results showed in Fig. [2] 
suggest slightly different effective temperatures along y 
and z directions. However, the difference between these 
temperatures is within error bars of our simulations. 

Derivation. — We make two assumptions about the 
external potential acting on the tagged particle. First, we 
assume that it is transverse, i.e. that it varies only in 
the plane perpendicular to the flow, d^ ■ TV'''^^{ri) = 0. 
Second, we assume that the potential is slowly varying. 
To make the latter assumption explicit and to facilitate 
the subsequent gradient expansion we introduce a small 
parameter e and write the external potential as V'^^^ (eri ) . 
As usual, at the end of the calculation we will put e = 1. 
We emphasize that no assumption is made regarding the 
strength of the external potential. 

It is natural to assume that for a slowly varying, trans- 
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Fig. 2: (Color online) Comparison of shear rate dependence 
of effective temperatures obtained from fits of tagged particle 
steady state density distributions to Eq. (|6]) and from the vio- 
lation of the Einstein relation, (a) velocity gradient direction, 
y. (b) vorticity direction, z. 

verse external potential the tagged particle density distri- 
bution n^i will also be slowly varying and transverse. To 
find a steady state equation for n^ we use the Chapman- 
Enskog [551 expansion. Specifically, we follow Titulaer [53] 
who applied the Chapman-Enskog expansion to the prob- 
lem of deriving the so-called contracted, Smoluchowski de- 
scription of a Brownian particle in position space from 
the full Fokker-Planck description in velocity and position 
space. We note that for linear equations considered in 
Ref. [53] and in our Letter, the Chanipan-Enskog proce- 
dure is closely related to a variant of perturbation theory 

Before applying the Chapman-Enskog expansion, we 
have to face the fact that, although the tagged particle 
density n^ is slowly varying, due to the inter-particle inter- 
actions TV-particle distribution P^ is not a slowly varying 
function of the tagged particle position ri. To disentan- 
gle the slow variation of P^ with ri , which is induced by 
the external potential, and the rapid variation P^ with ri, 
which originates from the inter-particle interactions, we in- 
troduce new variables: Ri — eri, R2 — r2i, ..., Rat = r^vi. 
In these new variables the steady state equation ([T]) has 
the following form: 



r!(0) + ef](i) -I- £2^(2) 
In Eq. ([n]). 



P3^(Rl,...,RAr)=0. (11) 



N 



N 



^^"' = -5r, ■ r ■ Ri -f ^ ^r, • Dod-Ei^ +DoJ2 ^^^ 

1=2 \ j=2 

N N \ 

-fio J2 F(R,,)+Mo5lF(-Rj)-r-R« (12) 

i^j=2 ]=2 J 

CN N \ 

i^o^^R. +Mo^F(-RO (13) 

i=2 i=2 / 

N N 

- 5] 9r. • MO (dn^V^'^'i-Ri)) - ^0 ^ ^r, • Or, 

i=2 i=2 

^('' = 9ri • {D^d^, + (9R,y'='='(Ri))) . (14) 



We now follow Titulaer [53] and simultaneously find a 
special solution of Eq. PT|) 

P,^(Rl,...,RAr) = ns(Ri)PiO)(R2,...,RAr) (15) 

-fePW(Ri,R2,...,Rw) + e'PfHRi,R2,...,Rw) + ... 

and a steady state equation for the tagged particle distri- 
bution, 



(p(0) + gp(l) + g2p(2) ^ J| ^^(j^^) ^ 0^ 



(16) 



Here Eq. (fTB|) is obtained from integrating Eq. ([TT|) over 
Ri, i> \. Thus, e.^., 

p(°)n,(Ri) = /"dR2...dRArr!(o)n,(Ri)Pf)(R2,...,Riv), 

(17) 

and 

p(i)n,(Ri) - /dR2...dR7vf^^''r^,(Rl)Pio)(R2,...,R7v) 

+ /"dR2...dRivfi^°'Pi'HRi,R2,...,Rjv). 

(18) 

It should be emphasized that in the Chapman-Enskog type 
procedure the steady state tagged particle distribution n^ 
is not expanded in e. Moreover, higher order terms Pg , 
i > 0, are determined up to an arbitrary multiple of Pg 
[53] . We follow the usual procedure [23j and impose con- 
ditions 

/"dR2...dRArPi')(Ri,...,RAr) == for i > 0, (19) 

which fix these contributions to be zero. It should be 
noted that due to these conditions the tagged particle dis- 
tribution rig is completely determined by the first term in 
expansion ([T5|) . 

We substitute p5|) into steady state equation pTJ) . solve 
iteratively for Ps , and calculate I?'-*^. The results of 
this procedure are as follows. First, Pg is the transla- 
tionally invariant solution of the Smoluchowski equation 
with shear flow in the absence of an external potential. 
This solution depends only on R^, i > 1, and is normal- 
ized in the following way: J rfR2...dR7vPg = 1. Second, 
25(0) — Q-^^ -F-Ri . Thus, for ng(Ri) that varies only in the 
plane perpendicular to the flow P^") does not contribute to 
steady state equation p6)) . Third, T)^^^ vanishes. Fourth, 
I?'-2) is determined by the following equation 

I?(2)„(R^) = (20) 

{D,,dR, ■ dn, + dn, ■ Mo {dn,V'^\Ri))] n(Ri) 

N 
-Or,- dR2...dRAr/io^F(-R,)Pi^HRl,-,Riv). 



i=2 
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To complete the calculation, 
terms of order e in Eq. (TTl]) . 



we determine P} using 



17(°)pii)(Ri,...,R^) = -r!(i)n3(Ri)Pi°)(R2,...,R 



N) 



where Fi is the force acting on the tagged particle, Fi = 
^j^2r(i"ij)- Following very similar steps one can show 
that 
(21) /x = M + M'" (27) 



Using the explicit form of il^^\ Eq. (fT3)) . it can be showed 
that the source at the right- hand- side of Eq. (f2T|) consists 
of two terms, which are proportional to [9R.jns(Ri)] and 
[9R,y-*(Ri)]«s(Ri): 

-17(i)n,(Ri)pW(R2,...,RAr) 

= [aR,n,(Ri)] 2i?o^5R, +^o^F(-R,) 

\ «=2 jjtl 

xPPCR2,...,'Rn) 

N 

+ [dR.V^^'iRi]] n,(Ri)Aio^aR,Pi°)(R2,...,Rjv). 

(22) 

To get Ps we now have to solve Eq. (121]) . While in- 
verting operator fl^°' we note that since both ns(Ri) 
and y"'*(Ri) vary only in the plane perpendicular to the 
flow, operator [f2^°)] will not act on them. Therefore, 

Ps also consists of two terms which are proportional to 
[5Rins(Ri)] and [aRiy''''*(Ri)] ^^(Ri). After these two 
terms are substituted into Eq. (^0)1 we can easily derive 
the following equality 

2?(')n(Ri) = [^R, • D • 5r, (23) 

+dR,-ti-{dR,V°^\Ri))]niR,). 

We now show that tensors D and /x, which enter into Ec^. 
(|23|) . are the self-diffusion tensor and the tagged particle 
mobility tensor, respectively. Explicitly, 



D = P>oI + D" 



(24) 



where I is the unit tensor and the interaction contribution 
D™ is given by the following formula, 

N 

D'" = / dR2...dRAr^F(-R,)[l](°)]^i 



N 



i=2 j^l 

cP,W(R2,...,Rjv). 



(25) 



After returning to the original variables ri,...,rjv and 
putting e = 1, the interaction contribution D™ can be 
written in the following form: 



D'" = ^ I dri...drNf^oFinj^i2Dodr,-fioFi) 



xPio)(ri,...,r^.) 



(26) 



where 

m'" -- ^ J dri...drNfioFin;^Dod,,P^°\ri,...,rN). (28) 

The above expressions are identical to the expressions for 
the self-diffusion tensor and the tagged particle mobility 
tensor that have been derived in Refs. [^[TT|. 

Finally, we return to the original variables in Eq. (P5|) . 
put e = 1, use the fact that both the external potential 
and the tagged particle density vary only in the plane per- 
pendicular to the flow, and arrive at Eq. ([SJ. 

Discussion. — We showed that the effective temper- 
ature defined through the violation of the Einstein rela- 
tion determines the tagged particle density distribution 
induced by a slowly varying external potential acting on 
the tagged particle. As a consequence, the same effective 
temperature enters into the violation of the Einstein re- 
lation and the modified static linear response relation for 
the tagged particle. It would be very interesting to ex- 
tend both results to at least some collective properties. 
This would provide additional fundamental justification 
for introducing the effective temperature. It would also 
allow some deeper theoretical understanding of the very 
intriguing observation made by Liu and collaborators [5] 
that for some combinations of variables one should use dy- 
namic linear response relations to determine the effective 
temperature whereas for other combinations of variables 
one could use static linear response relations. 

We note in this context that the most obvious (and per- 
haps the most naive) generalization of the present results 
to the collective density is not valid already in the linear 
response regime: we showed earlier [7] that the effective 
temperature obtained from the collective version of the 
modified static response relation does not agree with the 
effective temperature defined through the violation of the 
Einstein relation. 

It would also be very interesting to extend the present 
results, which pertain to non-linear response to an exter- 
nal potential, to material properties, in the spirit of the 
numerical investigation of Haxton and Liu |15| . 

The present results (and essentially all previous work on 
fluctuation-dissipation relation violation in systems of in- 
teracting particles) are concerned with a fluid undergoing 
a simple shear flow. A flow with a non-uniform velocity 
gradient, e.g. a Poisseuile flow, would introduce another 
length scale into the problem. Having three length scales 
in the problem, a microscopic length (colloidal particle 
size), the characteristic length of the external potential 
and a length characterizing the non-uniform flow, would 
make the present approach inapplicable. 

Finally, we should note that, while we derived a formal 
microscopic expression for the effective temperature, Eq. 
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O together with Eqs. (PHI^ and (PTll^ . we did not 
calculate the effective temperature theoretically. To do 
this, one would have to develop a theoretical approach to 
calculate interaction contributions to the self-diffusion ten- 
sor (|26p and the friction tensor (1^51) . In an earlier work [9 
one of us derived (but did not evaluate) mode-coupling ex- 
pressions for these quantities. A different mode-coupling 
approach was proposed by Kriiger and Fuchs [IT]. These 
authors evaluated their expression for the effective temper- 
ature only in the vicinity of the mode-coupling transition. 
This fact makes a direct comparison of our simulational 
results with their predictions difficult. 

We thank E. Flenner for comments on the manuscript 
and gratefully acknowledge the support of NSF Grant 
CHE 0909676. GS thanks Yukawa Institute, Kyoto Uni- 
versity, where this work was completed, for its hospitality. 
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